Magnetization dynamics induced by in-plane currents in ultrathin magnetic 

nanostructures 
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Ultrathin magnetic systems have properties qualitatively different from their thicker counterparts, 
implying that different physics governs their properties. We demonstrate that various such properties 
can be explained naturally by the Rashba spin-orbit coupling in ultrathin magnetic systems. This 
work will be valuable for the development of next generation spintronic devices based on ultrathin 
magnetic systems. 



Electric control of magnetic systems carries high po- 
tential towards device applications [l|,Q such as magnetic 
memory and logic. Spin-transfer torque (STT) is an 
efficient way to achieve the electric control of magnetic 
nanostructures. In view of device applications, magnetic 
nanostructures such MgO-based magnetic tunnel junc- 
tions are superior to silicon-based nanostructures in the 
simultaneous realization of nonvolatility and speed, but 
are estimated to require 100 times more energy Q than 
silicon-based CMOS devices to write an information bit. 
This energy cost problem limits scope of device appli- 
cations based on magnetic nanostructures. Since the 
writing energy decreases as a magnetic layer in magnetic 
nanostructures becomes thinner properties of ultra- 
thin magnetic layers are under intense investigation 

While the magnetization switching for the information 
writing is conventionally achieved by a current perpen- 
dicular to a magnetic layer, a recent experiment [7[ found 
that an in-plane current can also switch the uniform mag- 
netization of an ultrathin (ps 1 nm) magnetic layer (Co) 
sandwiched between a heavy metal layer (Pt) and an 
oxide layer (AlO^;) (Fig. [1]). Since the cross-sectional 
area (in yz-plane) for the in-plane current can be orders 
of magnitude smaller than the cross-sectional area (in 
xy-jAane) for the perpendicular current, this alternative 
switching scheme may reduce the current required for the 
switching and the switching energy. It was also reported 
that the magnetic domain wall (DW) in the ultrathin 
magnetic layer moves as fast as 400 m/s Q when in-plane 
current is supplied. This velocity is about 4 times higher 
than the highest velocity reported for thicker magnetic 
layers Thus the in-plane current effects on ultrathin 
magnetic systems open an attractive alternative path to- 
wards powerful spintronic devices. 

Ultrathin magnetic systems are interesting in view of 
fundamental science as well. Various features of ul- 
trathin magnetic systems cannot be explained by ex- 
isting theoretical knowledge learned from measurements 
on thicker counterparts; for example (i) magnetization 
switching by in-plane current 0] instead of perpendic- 



FIG. 1; (color online) Schematic structure of ultrathin mag- 
netic nanostructure, where an ultrathin {~ 1 nm) ferromag- 
netic layer is sandwiched between a heavy metallic layer and 
an insulating oxide layer. Examples include Pt/Co/AlOi and 
Ta/CoFeB/MgO. 



ular current, (ii) DW motion against the electron flow 
direction instead of along it, and (iii) anomalously 

high DW speeds Q. These anomalies imply that ultra- 
thin magnetic systems are not mere thin limits of thicker 
counterparts but are qualitatively different systems gov- 
erned by different physics. A clear understanding of their 
core physics will be highly valuable for developments of 
spintronic devices. 

Since the upper and lower layers of the ultrathin mag- 
netic layer are made of quite different materials (Fig. [1} , 
the ultrathin magnetic layer has the structural asymme- 
try. Here we demonstrate theoretically that if the Rashba 
spin-orbit coupling (RSOC) [llj due to the asymmetry 
is sufficiently strong, all the anomalies (i), (ii), and (iii) 
can be explained naturally. To be specific, we show that 
the in-plane current gives rise to a torque proportional 
to Q!Rm X [m X (z X j^)] [l^, where q;r is the param- 
eter describing the strength of the RSOC, jg is the in- 
planc current density in the ultrathin magnetic layer, z 
is the unit vector perpendicular to the layer, and m is 
the unit vector along the magnetization in the ultrathin 
magnetic layer. Since this torque has the same form as 
the Slonczewski STT Q in magnetic multilayers, we call 
it Slonczewski-like STT (SL-STT). We demonstrate that 
the SL-STT explains all three anomalies naturally. Re- 
cently an experiment [l3| proposed that the anomaly (i) 
may arise due to the spin Hall effect (SHE) in the heavy 
metal layer. Here we demonstrate that the SHE alone 
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cannot explain the anomalies (ii) and (iii) . 

First, we demonstrate that RSOC generates the SL- 
STT. Conduction electrons in the ultrathin ferromag- 
netic layer can be described by the Hamiltonian, 
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where p^p is the momentum operator, rUe is the elec- 
tron mass, Jex (> 0) is the exchange energy, and cr 
is the Pauli matrix. The third term in Eq. ([T]) is the 
RSOC Hamiltonian and the last term describes spin re- 
laxation processes such as electron scattering. When 
expressed in terms of the kinematic velocity operator 
Vop = Pop/me + (q;r/?i)z X <T, Eq. ([T]) becomes H = 
'^eVop/2 -|- Joxf ■ m + H^q\, where a trivial c-number has 
been neglected. Due to the strong exchange energy Jex, 
which is the largest energy scale affecting the spin dynam- 
ics in conventional metallic ferromagnets (such as Co or 
Fe), electrons in energy bands still carry largely major- 
ity and minority spin directions. However within each of 
these bands, the individual spins and the net spin are not 
nessarilily coUinear with the magnetization due to tilting 
of the spins by RSOC. In the following, we treat the elec- 
trons on the majority and minority sheets of the Fermi 
surface separately, but need to treat each spin density as 
a vector s± = {<y)±, since it is not aligned with the mag- 
netization. Here (• • • )± denotes local average over ma- 
jority/minority electrons. The spin continuity equation 
determined by H allows one to derive the Bloch equation 
for s±. 
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where Tex = fif^Jex, v± = (vop)±, and T = [cr, i^rcil/ift, 
and the spin-current tensor density J^± ~ v± (g) s± . The 
approximation ((<T)»(voo),-|-(vop),-(cr),)+ = 2(s±)j(v±)j 
is used in Eq. ^ RecaUing that STT is deter- 

mined by the transverse component 5s± of s-t perpen- 
dicular to m, it is useful to separate s± into longitudinal 
and transverse components, s± — =pn±m + 5s±, where 
n± = =Fs± ■ m is the longitudinal spin accumulation. 
One then makes the relaxation time approximation [l5| , 
(r}j- = — (Js-t/Tsj., where t^^ is the transverse spin relax- 
ation time. In this approximation, the relaxation of the 
longitudinal spin component is neglected since the trans- 
verse relaxation is much faster in conventional metallic 
ferromagnets and also the longitudinal spin component 
does not affect the STT. When v± is assumed to be ho- 
mogeneous [lij within the magnetic layer, one can deter- 
mine from Eq. ([2]) the following approximate solution for 
(5s±, 



(5s± = ±n±Tc 
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where /? = t^x/tsi and Df = 9t -I- v± • V. Correc- 
tions to this solution are of higher order in t^x and 
may be neglected in the strong Jgx limit. Finally from 
the relation T = ^bt^x^ ^ between STT T and 
the total transverse spin density = 5s+ + (5s_, one 
obtains the Landau-Lifshitz- Gilbert equation dlsA/dt = 
-70M X Heff + {ao/Ms)M X dM/dt + T, where Heg is a 
sum of an external magnetic field and effective magnetic 
fields due to magnetic anisotropy and magnetic exchange 
energy. M = Msin is the magnetization and Ms is sat- 
uration magnetization in the ultrathin magnetic layer. 
After grouping together terms of the same structure, one 
obtains 
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where Hr is the additional effective magnetic field due to 



RSOC, jg = — e(n+v+ + ?7._v_), and Pjg 



-e(n_| 



n_v_). Note that the first and fourth terms on the 
right-hand-side of Eq. ([4]) contain the renormalized gy- 
romagnetic ratio 7 and the renormalized Gilbert damp- 
ing a given by 70/7 = 1 + (71+ - rt_)/[Ms(l -I- /32)], and 
7q;/7o = Q;o + /?(n+-n_)/[A/s(l + /3^)]. The last and sec- 
ond to the last terms are the nonadiabatic STT [l^ [l3| 
and the adiabatic STT [l^ • These four terms govern the 
magnetization dynamics in thicker magnetic systems. 

The second and third terms are additional STTs due 
to RSOC and may have sizable magnitude only in ultra- 
thin magnetic systems since q;r decays as the magnetic 
layer becomes thicker. The second term — 7M x Hr has 
the same structure as the first term and thus we call it 
field-like STT (FL-STT). The FL-STT was derived pre- 



viously [19|, l20| and Ref. l21| reported its experimental 
confirmation. The FL-STT alone, however, can explain 
none of the anomalies, as we demonstrated recently [2^ . 
The third term is the very SL-STT that this calculation 
aims to derive. 

Next we explain the anomalies in terms of the SL-STT. 
Explanation for the anomaly (i) is trivial. Recalling that 
the SL-STT due to the in-plane current has the exactly 
same structure as the Slonczewski STT Q due to the per- 
pendicular current in metallic spin valve systems and that 
the Slonczewski STT induces the magnetization switch- 
ing in the spin valve systems, it is easy to understand 
that the in-plane current can induce the magnetization 
switching through the action of the SL-STT. 

The anomalies (ii) and (iii) are less trivial to explain 
since they arise from the combined action of the SL- 
and FL-STTs. To explain the anomalies, we begin with 
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FIG. 2: (color online) Two possible structures of the Bloch 
DW when the current flows to the right and a > /3 [2^ . Thick 
and thin arrows represent the directions of m and STTs, re- 
spectively. Note that the direction of m deviates from ±y, 
which is a generic feature of a moving DW [IB, [13 with a ^ 13. 



two possible structures of the Bloch DW (Fig. [5]) in ul- 
trathin magnetic systems (such as Pt/Co/AlOx [1, HH) 
with the perpendicular magnetic anisotropy. In the ab- 
sence of RSOC and when jp is not large enough to cause 
the DW structural instability (Walker breakdown) [l7j . 
the two structures are equivalent in term of both sta- 
bility and dynamics (same DW velocity wdw)- The 
first effect of RSOC is to break the dynamic equiva- 
lence; at the DW center in Fig. M^a)/{h), the SL-STT 
is anti-parallel/parallcl to the nonadiabatic STT, ef- 
fectively cancelling/enlarging the effect of the nonadia- 
batic STT. Recalling that the nonadiabatic STT deter- 
mines udw [3 [13 1 this implies that the SL-STT re- 
duces/increases vuw- When RSOC is sufficiently strong, 
it is even possible that the SL-STT over-cancels the nona- 
diabatic STT in Fig. Elja), so that wdw reverses its sign. 
When RSOC is even stronger, udw will be large with the 
reversed sign, implying that the DW in Fig. ^a) moves 
fast against the electron flow direction. By the way, the 
FL-STT does not affect wdw [llj since it is perpendic- 
ular to the nonadiabatic STT at the DW center. The 
second effect of RSOC is to break the stability equiv- 
alence. The effect of the FL-STT on the stability can 
be understood from the direction of the effective field 
Hr. Since the effective energy density — Hr ■ M is neg- 
ative/positive at the DW center for the DW structure in 
Fig. [2i;a)/(b), the FL-STT makes the DW structure in 
Fig. ^a.) more stable than the other structure. More- 
over when Hpt is sufficiently strong, the DW structure 
in Fig. [2jb) becomes unstable and evolves to the stable 
DW structure in Fig.[2i;a) By the way the SL-STT 
has much weaker effect on the stability than the FL-STT 
since, according to Eq. (|4]), the SL-STT is smaller than 
the FL-STT in magnitude by factor (3 and this nonadi- 
abaticity parameter /3 is known [Tij to be smaller than 



1. Then combining the above information, we find that 
there is only one stable DW structure [Fig. [2][a)] when 
RSOC is sufficiently strong and that it moves fast against 
the electron flow direciton, explaining both the anomalies 
(ii) and (iii). 

The experiment 0, which reported the anomalies (ii) 
and (iii), provides sufficient information to test this ex- 
planation. According to Eq. the magnitude of the 
SL-STT is 7/3|MxHr| = 7/3Ms|Hr| sin0, where (/) is the 
DW tilting angle between +y and M at the DW center, 
and the experiment reported |Hr| ^ IT at |jg| = 10^^ 
A/m^. On the other hand, the magnitude of the nona- 
diabatic STT is l3bj\dM/dx\ - (3b.jM^/\ where the DW 
width A is about 5 nm and hj = /iB-Pje/[eAfs(I -I- 0^)] 
is of the order of 100 m/s at |jj = lO^^ A/m^. Thus 
the relative magnitude of the SL-STT with respect to 
the nonadiabatic STT is roughly given by ctR sin 0, 
where the dimensionless parameter Q!r = TrajimeX/h^ = 
(7r/2)(7;3Ms|HR|)/(/3&jMs/A) is of the order of 10. Thus 
unless (p is too smafl, the SL-STT can be indeed larger 
than the nonadiabatic STT and reverse the sign of ?;dw- 

Figures ^a) and (b) show micromagnetic simulation 
results of Eq. ([4]) for the stable DW structure in Fig.[2][a). 
The inset in Fig. ^a) shows wdw as a function of q;r at 
fixed je = X • jg = +3 X 10^^ A/m^ (amounting to bj = 
+25 m/s). Note that as or increases, vnw changes its 
sign from negative (along the electron flow direction) to 
positive (against the electron flow direction). The main 
panel in Fig. ^a) shows ?;dw as a function of je at two 
fixed values of Sr, (black squares) and 10 (red circles). 
For QfR = 10, vuw changes from negative to positive 
at je ~ 3.5 X 10^^ A/m^ and goes above -1-500 m/s for 
je > 1.5 X 10^^ A/m^. Thus both the anomalies (n) 
and (iii) can be explained by RSOC if Q!r is sufficiently 
larger than I. By the way, results for the DW structure 
in Fig. [2jb) are not shown since, when Sr = 10, it is 
unstable and switches to the DW structure in Fig. [H^a) 
for je > 7.4 X 10^° A/m2. 

Next we demonstrate that the above explanation for 
the anomalies (ii) and (iii) does not work if the FL-STT 
is absent or very small. A recent experiment [isj on a 
somewhat different ultrathin magnetic system Ft /Co/ Al 
reported that the SL-STT has a sizable magnitude but 
the FL-STT is negligibly small. This situation has been 
attributed [l^ to a perpendicular spin current in Co gen- 
erated from an in-plane charge current in the heavy metal 
layer (Ft) through the SHE. In the absence of the FL- 
STT. the stabilities of the two DW configurations in 
Fig. [2] are governed by the SL-STT. Recalling that the 
DW anisotropy energy of the Bloch DW favors (p ~ 
and ±180°, and disfavors ±90°, it is evident that the 
SL-STT tends to destabilize/stabilize the DW structure 
in Fig. (la)/ (b). Thus when the SHE-induced SL-STT 
is sufficiently strong, the DW structure in Fig. EJb) is 
the only stable structure, which is opposite to the choice 
made by the FL-STT. Micromagnetic simulation for this 
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FIG. 3: (color online) Micromagnetic simulation results of the 
current-driven DW motion. In (a) and (b), RSOC effects are 
examined by using Eq. (j4| . vow (a) and the DW tilting angle 
(j> (b) [measured clockwise from +y direction in Fig. [2ja)] as a 
function of je for ur = (black squares) and 10 (red circles). 
As shown in (b), the Walker breakdown is suppressed when 
i3r — 10 whereas it occurs for > 1.0 x 10^^ A/m^ when 
Or = 0. The inset in (a) shows i^dw as a function of an for 
je = 3.0 X lO^'^ A/m^. All results in (a) and (b) are for the sta- 
ble DW configuration in Fig.[2Ka). In (c) and (d), SHE effects 
are examined by setting Hr = and instead adding to Eq. (|4} 
the SHE-induced SL-STT 7M x [(6»shM/Ms) x (Hsy)], where 
^^s = hje,N / {2eAdstF), tp is the thickness of the ultrathin 
magnetic layer, and je.N is the charge current density in the 
heavy metal layer. je,N = je is assumed, vbw (c) and (j) (d) as 
a function of je for ^sh ~ (black squares) and 0.048 [red cir- 
cles and green diamonds for the DW structures in Figs. [2]^ a) 
and (b)]. Red circles are not visible for je > 5 x 10^^ A/m^ 
since they overlap with green diamonds. The inset in (c) 
shows vuw as a function of ^sh at — 3.0 x 10^^ A/m^ for the 
stable DW configuration in Fi g. [2f b). The parameters for the 
simulation; a = 0.5, /3 = 0.25 Ms = 5.0 x 10^ A/m, the 
perpendicular magnetic anisotropy constant Ku ~ 1.0 x 10^ 
J/m^, the exchange stiffness constant Aa — 1.0 x 10^^^ J/m, 
P = 0.7, 7/27r = 28.0113 GHz-T~\ and tp = 0.6 nm. 



stable DW structure indicates [inset in Fig. ^c)] that 
vdw does not change its sign as 0sh grows, which is in 
contrast to the RSOC effect [inset in Fig. (Ha)]. Here 
^SH is the spin Hall angle representing the strength of 
SHE. This result is natural since for the stable DW struc- 
ture, the SL-STT is parallel to the nonadibatic STT at 
the DW center. Note also that large enhancement of 
|wdw| does not occur [Fig.tSjc)] either, since the SL-STT 
in Fig. E^b) suppresses the deviation of (f> from ±180°, 
thereby suppressing its own magnitude, which is pro- 
portional to sin(/). Thus SHE alone cannot explain the 
anomalies (ii) and (iii). By the way, the DW structure 



in Fig. mja) may exhibit the sign-reversed [red cir- 
cles in Fig. [3^c) near — 4: x 10^^ A/m^] right before 
it loses its stability at je ~ 5 x 10^^ A/m^ [downward 
jump of red circles by 180° in Fig. EJd)]. However this 
sign-reversal should be contrasted with the sign-reversal 
observed Q in the current range, where only one stable 
DW structure exists. 

To conclude, we presented the theory of the RSOC- 
induced STTs, which explains various anomalous fea- 
tures of the magnetization dynamics in ultrathin mag- 
netic systems. It will be valuable for future works 
to utilize ultrathin magnetic systems for next genera- 
tion spintronic devices. Important future research di- 
rections include exploring other material combinations 
such as Ta/CoFeB/MgO, utilizing both RSOC and SHE 
to achieve enhanced properties, and understanding the 
dependence of the RSOC strength on magnetic layer 
thickness and material combinations. It is also inter- 
esting to explore possible connections between this work 
and other interesting phenomena in ultrathin magnetic 
systems such as magnetization control via electric volt- 
age [2^ and strong perpendicular magnetic anisotropy Q 
in Ta/CoFeB/MgO systems. Lastly, during the prepara- 
tion of this manuscript, we received a calculation [25[, 
which also presents a derivation of the SL-STT but does 
not discuss implications of RSOC on the DW motion. 
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